An overlapping Schwarz method for spectral element simulation of three-dimensional incompressible ows

As the characteristic propagation speed is innnite for unsteady incompressible ows, solving for the pressure operator is the most computationally challenging phase of ow simulation. EEcient solution is dependent upon the availability of fast solvers for sparse linear systems. In this paper we develop an Schwarz preconditioner using overlapping subdomains for the pressure. These local subdomain problems are derived from tensor products of one-dimensional nite element discretizations and admit use of fast diago-nalization methods based upon matrix-matrix products. In addition, we use a coarse grid projection operator applied directly to the pressure on the interior Gauss points whose solution is computed via a fast parallel direct solver. The combination of overlapping Schwarz preconditioning and coarse grid solver provides as much as a four-fold reduction in simulation time over previously employed methods based upon deeation for parallel solution of multi-million grid point ow problems.